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STELLAR POPULATION GRADIENTS IN ULIRGS: IMPLICATIONS FOR GAS INFLOW TIMESCALES 
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ABSTRACT 

Using longslit, optical spectra of Ultraluminous Infrared Galaxies (ULIRGs), we measure the evo- 
lution in the star-formation intensity during galactic mergers. In individual galaxies, we resolve kpc 
scales allowing comparison of the nucleus, inner disk, and outer disk. We find that the strength of 
the H/3 absorption line increases with the projected distance from the center of the merger, typically 
reaching about 9 Aaround 10 kpc. At these radii, the star formation intensity must have rapidly de- 
creased about 300-400 Myr ago; only stellar populations deficient in stars more massive than Type A 
produce such strong Balmer absorption. In contrast, we find the star formation history in the central 
kpc consistent with continuous star formation. Our measurements indicate that gas depletion occurs 
from the outer disk inwards during major mergers. This result is consistent with merger-induced gas 
inflow and empirically constrains the gas inflow timescale. Numerical simulations accurately calculate 
the total amount of infalling gas but often assume the timescale for infall. These new measurements 
are therefore central to modeling merger-induced star formation and AGN activity. 
Subject headings: galaxies: starburst — galaxies: evolution — galaxies: active — galaxies: formation 



1. INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs) are some 
of the most luminous objects in the local universe with 
log(L/i^/LQ) > 12 ~ a result of starburst and nuclear ac- 
tivity triggered by major mergers. The number density of 
ULIRGs per unit comoving volume has evolved strongly 
causing ULIRGs to become a rare phenomenon locally, 
but infrared-luminous sources appear t o dominate the 
star-f orming activity beyond z ~ 0.7 (|Le Floc'h et al.l 
120051 ). By studying local analogues of these important 
objects we can better understand the physical processes 
regulating their evolution. 

When stellar disks merge, phase-space density remains 
constant and therefore gen erally too low to m atch that 
of an elliptical galaxy (.Hernguist et al.l I1993D . Merg- 
ing gas-rich galaxies, however, allows dissipation; and 
the starburst increases the central phase-space density 
([Robertson et al.l 12006) . The extremely high star for- 
mation rates pre sent in ULIRGs i ndicate high central 
gas densities (Sa nders et al.l I1986D . so ULIRGS meet 
the conditions for dissipative mergers. Further evidence 
that ULIRGs mark this morphological transition comes 
from their stellar kinematics and surface brightness pro- 
files, which resemb le those of field elliptical galaxies 
([Genzel et al.ll200lD . In the ULIRG phase, however, the 
gas mass is much higher than that of an elliptical galaxy. 
The gas fraction blown away by stellar feedback is un- 
clear, but half or more is likely consumed by star forma- 
tion. Where in the galaxy this star formation occurs is 
relevant to the steepening of the stellar age ~ metallic- 
ity gradient and the fueling of the central supermassive 
black hole in the merger remnant. 

The high central concentration of molecular gas in 
ULIRGs indicates a strong bu rst of star formation in the 
central region of a kpc or less ([Sanders et al.llT986[ ) . The 
rotation measured in some outflows, however, suggests 
a larger region drives a galactic wind (Martin 2006). In 
this paper, we use longslit, optical spectroscopy to exam- 



ine radial gradients in the stellar populations of ULIRGs. 
Integrated spectra, in contrast, are heavily weighted to- 
wards the recent star formation activity in the central 
kiloparsec. We present the sample in Section 2. In Sec- 
tion 3, we measure radial gradients in stellar spectral di- 
agnostics. We focus on the strength of Balmer absorption 
and emission. The hydrogen recombination rate directly 
measures the ionization rate and is sensitive to the num- 
ber of high-mass (and therefore young) stars. Balmer 
absorption is strongest in Type A stars, and high val- 
ues indicate a paucity of OB stars and declining star 
formation rate. In Section 4, we associate the cessation 
of star formation with depletion of the local gas den- 
sity. We discuss how our result constrains the gas inflow 
timescale and compare our measurement to the dynam- 
ical timescale. Section 5 summarizes our conclusions. 

2. DATA 

2.1. Observations 

The Keck II echelle te spectroscopy of 2 J y ULIRGs, de- 



r oscopy 01 z J y 
scribed previously in iMartinI ([2005L [2006) , is well suited 
to our study of H Balmer emission/absorption. These 
data resolve scales of 0.8", about 1.5 kpc, spatially and 
~ 70 km s~^ spectrally. The 20" slit is long enough to 
sufiiciently sample the sky at both ends of the slit, en- 
suring accurate sky subtraction and continuum intensity. 
The slit position angle includes a second nucleus when 
present (Table [T]). We inspected 43 ULIRG spectra and 
selected 25 for analysis of any stellar II/3 absorption. This 
subsample presents II/3 emission such that the full width 
of the emission line in their nuclear spectrum, defined by 
the wavelengths where the trough changes from emission 
to absorption, is no more than 1800 km s^^. 

The resulting subsamp le includes a r ange of stages in 
the merger progression. iVeilleux et all (|200g ) classified 
the merger s tage for 11 of the 25 o bjects. Using K-band 
images from iMurphv et al.l ([19961) , and their identifica- 
tion as multiple nuclei systems, we apply the Veilleux 
scheme to the remainder finding 8 pre-mergers (i.e. two 
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band morphology. "S", "D", and "M" denote a sin- 
gle, double and mult iple nucleus morphology respectively 
(Murphy et al.|[T996l ). Col. (3): Optical R band morphol- 
ogy. Similar to R band morphology with an additional 
"t" to signify an extended tidal region, and "W" to de- 



note a widely separated pair llMurphv et al.lfT996l). C ol. (4) 



Merger Classification based on lVeilleux et al ■l (|2002l l. This 
classification separates wide binary pre-mergers (Ilia), 
close binary pre-mergers (Illb), diffuse merger(IVa), com- 
pact merger (I Vb), and old rnergers (V). Col. (5): IR lu- 
minosity from (|Murphv et all 119961 '). Col.(6): redshift. 
Col. (7): Position angle of slit 

Classifications directly from I Veilleux et al.l (|2002f ) 

nuclei), a multiple merger (IRAS 10565+2448), 12 merg- 
ers, and 4 old mergers. The ULIRGS in the merger stage 
tend to present higher surface brightness, so our analy- 
sis of spatial gradients in Balmer absorption is heavily 
weighted towards this merger phase. 

2.2. Measurements 

We examined the positional dependence of the H/3 line 
profile using a sliding aperture for spectral extraction. 
We matched the aperture width to a seeing element 
(~ 0".8) and extracted a spectrum at each line of the 
spectrogram. At some spatial positions, as illustrated 
in Fig. [TJ the spectra present a two-component H/3 line 
profile. A broad absorption trough underlies the narrow 
emission line. We fit the absorption and emission com- 
ponents simultaneously using non-linear least-squares fit- 
ting with MPFIT in IDL. These spectra exhibit continu- 
ous change from position to position due to atmospheric 
smearing, so parameters found from fitting one aperture 
position were supplied as the guess for the fit of the next 
position. The small, gray points in Figure [2] show the re- 
sult. Over scales of a few kpc, this analysis often reveals 



an increase in H/3 absorption equivalent width, W^f/3, 
with distance from the peak continuum emission in the 
R band. 
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Fig. 1. — The aperture selected at 2 kpc north west from the 
center in IRAS17208-0014 is shown above with a fit to both Hp 
components simultaneously. Emission component is a narrow line 
contributed by photoionization of the ISM. The absorption compo- 
nent is provided by the presence of A stars in the stellar population. 
The emission line is narrow enough that the absorption line is not 
affected in the fit. 



To eliminate effects due to correlated errors, we ex- 
tracted spatially separated spectra based on the shape of 
these radial H/3 absorption and continuum surface bright- 
ness trends. These spectra are shown in Figure [31 where 
the larger black points in Figure [2] show the H/3 absorp- 
tion equivalent width obtained by fitting these spectra. 
We measured Ha surface brightness in these same aper- 
tures. The H(5 lines could only be fitted in more than 3 
apertures for one object. 

To correct the Ha emission for underlying stellar 
Ha absorption, we measured the H/3 absorption equiv- 
alent width and assumed the Ha absorption equiva- 
lent width was two-thirds as large. Inspection of syn- 
thet ic spectra computed for stellar population mod- 
els (jBruzual &: CharlotI 120031 ) verified that this scal- 
ing works over a broad range of star formation mod- 
els, including those with truncated star formation his- 
tories. The maximum value of the correction, 3A, is 
small compared to the emission equivalent width. The 
flux in Ha emission is also corrected for extinction. 
We determined the color excess, E(B - V), per aper- 
ture by measuring the Balmer decrement and using 
a standard reddening law (Ry — 4.05 — Av /E(B - 
V)) and the extinction law from iCalzetti et all (|2000). 
When significant, we discuss the correction for at- 
tenuation by Galactic dust on a galaxy-by-galaxy ba- 
sis in Section 13.21 We estimate areal star forma- 
tion rate using the conversion S'Fi?(M0 yr^^ kpc^^) — 
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Fig. 2. — key : Upper: The equivalent width of H/3 due to stellar 
absorption is plotted versus position. The light grey points are 
the measurements from the preliminary sliding aperture analysis 
of the spectrum. Black points are chosen apertures based on trends 
in the preliminary analysis. The solid black line shows continuum 
intensity near Ha as a function of position in units normalized 
to the peak of the continuum, illustrating the alignment with the 
images in the lower panel and an indication of the relative signal 
to noise in the extended regions. The solid grey line is the same 
for the continuum near H/3. Lower: An R band image (courtesy of 
T. Murphy) is shown with the position of the slit plotted to show 
where the 1" X 20" slit was placed. We can see increasing Y^up 
in the extended regions with a sharp drop in the nucleus. Figures 
2. a — 2.j are available in the online version of the Journal. 

[i(iJa)/A(kpc2)]/1.26 X lO-^i erg s"! (|Kennicutdll998l) . 
3. RESULTS 

The presence of H/3 emission and absorption reveal dif- 
ferent epochs of star formation. Recombination emission 
traces ionization by very massive stars, and Balmer ab- 
sorption is strongest in Type A stars. Strong Balmer 
absorption dominates the stellar spectrum when hotter 
stars are absent. Figure S] shows how the measured H/3 
absorption equivalent width changes along the slit. The 
projected distance is measured from the position along 
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Fig. 3. — The extracted spectra used for measurement are plotted 
in grey with the fit to the absorption overplotted in black. Each 
of these spectra correspond to one of the black points in Figure [2] 
Figures 3. a — 3.j are included in the online version. 



the slit with the highest R band surface brightness. This 
center typically coincides with the position of the nucleus 
measured in the K band. IRAS17208-0014 has two max- 
ima in R band continuum, however only the brightest of 
these maxima coincides with a K-band nucleus, which 
we regard as the center of the object. IRAS 15245+1019 
has two nuclei resolved in K-band images (|Murphv et alj 
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Il996f) separated by - 4", but we use the nucleus with 
the highest R-band surface brightness to define the cen- 
ter. K-band imaging of IRAS20046-0623 does not resolve 
two nuclei, but has elongated structure over 4" from east 
to west. The position of the maximum continuum sur- 
face brightness near Ha differs from that of H/3as well. 
For this object we chose to retain the definition of center 
based on the maximum R-band surface brightness since 
the K-band centroid is not directly measured in these 
data. 

The strength of the stellar Balmer absorption typically 
increases with distance from the center. For example, in 
FigurelHthe H/3 absorption strength across IRAS 09111- 
1007 increases steadily out to 7 kpc. This behavior char- 
acterizes merger systems presenting relatively symmet- 
ric structure on either side of the nucleus. Objects such 
as IRAS20046-0623 observed at an earlier stage in the 
merger progression often present more complicated mor- 
phology. The spatial gradient in H/3 equivalent width is 
not smooth like that in IRAS09111-1007, but the absorp- 
tion strength is still largest away from the center. The 
surface brightness in the post mergers is generally too low 
to look for the effect at large radii. The galaxies where 
such gradients can be measured fall in or between the 
close binary pre-merger stage (HIa) to compact merger 
stage (IVb). The 10 targets caught at these stages form 
our spatially resolved subsample, shown in Figure |4] and 
discussed individually in Section 3.4. 
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Fig. 4. — Equivalent width of H/3 absorption is plotted as a func- 
tion of distance from the maximum continuum surface brightness. 
The absorption strength generally increases away from the center. 



Emission effects the line profile and equivalent width 
measurement less in Kd. However, IRAS 23365-f 3604 is 



the only galaxy where our data quality is sufficient to 
study their line profiles in more than 3 apertures. We fit 
the H/3 and US emission/absorption profiles simultane- 
ously, requiring the absorption features to have the same 
width and Doppler shift but allowing different normaliza- 
tions. We apply the same requirement to the fit of the 
emission components. Figure [5] shows that the equiva- 
lent widths of these lines increases with distance from the 
center. Hence our results are not strongly influenced by 
our decomposition of the emission and absorption com- 
ponents of the H/3 profile. We can use the measured 
strength of Balmer absorption to constrain the recent 
star formation history at different radii. The measure- 
ments in Figure [5] also show that the equivalent width 
of H/3 and US are similar in the apertures that could 
have both balmer lines measured. The similarity is in- 
consistent with the population modeling which predicts 
slightly greater equivalent widths for US. Dominance 
of spectral types AO through A5 could possibly explain 
this. A star formation history that allowed these stars 
to dominate might be required if measurements of more 
ULIRGs show this same inconsistency. 
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Fig. 5. — Equivalent widths of H/3and H<5 absorption are plot- 
ted as a function of distance from the central continuum peak in 
the chosen apertures for IRAS23365-I-3604. The Balmer equivalent 
widths appear to follow the same pattern. 



3.1. Truncation Timescales from Stellar Population 
Synthesis 

Using GAL AXE V (IBruzual fc CharlotI 120031) . we syn- 
thesized model spectra of the stellar population for three 
different star formation histories: (1) S function, (2) con- 
tinuous star formation history that is extinguished after 
40 Myr, and (3) continuing star formation (CSF). We 
assumed solar metallicities. We measured Wnfi in these 
spectra using the same techniques applied to the obser- 
vations. We illustrate the evolution of the Wh^ as the 
stellar population ages in Figure El Around 850 Myr, 
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continuous star formation reaches a maximum of 6 A 
in equivalent width as the population reaches an equilib- 
rium of births and deaths of A stars. In the other models, 
as O and B type stars die off, the stellar population dom- 
inating the stellar continuum emission gradually shifts to 
later stellar types. A single A5V star has W//^ near 20 
A, but the spread in stellar types does not allow the com- 
posite models to reach this limit. Once star formation 
is stopped, the W//^ increases steadily towards a maxi- 
mum of 9A at 400 Myr. After 400 Myr, the Wnp falls as 
the number of A stars declines, reaching 6A after about 
850 Myr. An equivalent width greater than 6 A clearly 
requires a very rapid decrease in the star formation rate 
in the past Gyr. 

The large values of found at radii R > 5 kpc 

in IRAS 23365-1-3604 are inconsistent with the contin- 
uous star formation rate model. Applying the relation 
between Whi3 and time since truncation, i.e. Figure |6] 
solid curve, we conclude that star formation ceased at 
least 100 Myr ago at i? > 5 kpc. The HS equivalent 
width varies much like H/3 in Figure |6] owing to its sim- 
ilar origin. The gradient continues towards the 
center of IRAS 23365-f 3604, and we interpret this as ev- 
idence that the time since truncation decreases towards 
the center. In other words, star formation is apparently 
being shutdown from the outside inwards. 

Another popular diagnostic of stellar age is the spec- 
tral break known as D4000. Our spectral bandpass does 
not completely cover this feature, but we can use the 
stellar population models along with our observed trend 
in Balmer equivalent width to predict the radial varia- 
tion in D4000. The D4000 index grows as the stellar 
population ages. Measurements of D4000 break the age 
degeneracy to either side of the maximum in Wh/s at 
400 Myr. Hence, we expect D4000 would show a steady 
increase with increasing radius in ULIRGs. 

Our simple stellar populations represent the most ex- 
treme limiting cases. A declining star formation rate 
produces a rising Wh/S similar to our models with com- 
plete cessation. The stronger the rate of change in the 
star formation rate, the closer the Whp comes to that of 
the truncated model. We can distinguish these two star 
formation histories by the absence/presence of Balmer 
emission in our ULIRG spectra. The presence of nebular 
emission indicates the level of ongoing star formation. 
Complete truncation of star formation activity reduces 
the ionization rate, and Balmer emission ceases. We ex- 
amine the radial depedence of the implied star formation 
rate galaxy by galaxy. 

3.2. Spatially Resolved Star Formation Histories 

In Section 3, we identified 10 galaxies having extracted 
spectra in a least four apertures with SNR > 5. Pre- 
merger and merger systems dominate this spatially re- 
solved subsample. We discuss the recent star formation 
history in each individually here. For reference. Table [2] 
lists the measured star formation rate and stellar popula- 
tion ages for different apertures in each object. The star 
formation rates near the center should be treated cau- 
tiously due to two important systematic errors. First, 
an active galactic nucleus may contribute to the ioniz- 
ing photon luminosity, thereby lowering the inferred star 
formation rate. Second, the reddening measured by the 



Balmer decrement in the central aperature reflects only 
the least obsured regions; the extinction may be too high 
in some regions for Balmer photons to escape causing us 
to underestimate the central star formation rate. These 
factors have little impact on our analysis of spectra ex- 
tracted outside of the central kiloparsec and our discus- 
sion of radial gradients. 

Our results suggest the following star formation histo- 
ries. 

IRAS15245+1019: This ULIRG has two K band nu- 
clei separated by « 4kpc. These two nuclei host star 
formation of 4.83 ±0.05 and 6.8 ±2.1 Mq yr'^ kpc^^ for 
the brighter south-east nucleus and dimmer north west 
nucleus respectively. The relative dimness in the visible 
band of the southeast nucleus is due to high levels of ex- 
tinction by dust. The brighter nucleus in the interaction 
has Wif^ 5.0 ± 0.1 A which indicates that star forma- 
tion has been going on for 0.22 ± 0.02 Gyr. At 7 kpc, 
however, there is little to no star formation occurring 
currently. W hp indicates that 0.2 ±0.1 Gyr has passed 
since truncation. 

IRAS20046-0623: From detailed analyses of galactic 
relative velocities and tidal tail characteristics, it is deter- 
mined that IRAS20046-0623 has g one through pericen - 
tcr in the last few times 10 years (jMurphv et al.l[200lh . 
The separated nuclei are masked within the extended 
central bulge that appears in both K band and R band 
imaging. The morphology of the extended tails are rep- 
resentative of two nearly orthogonal disks in the process 
of merging. Peaks in star formation (3.95 ± 0.05 and 
8. 3 ± 0.7 M(7) yr" " '" kpc~ ^) occur in the nuclei indicated 
in iMurphv et al.l ()2001[ ). As shown in figure [He, these 
nuclear regions have W^^^ that indicate star formation 
has continued for 0.18 ± 0.03 Gyr and 0.36 ± 0.07 Gyr. 
In the outer regions ~ 2 kpc beyond both nuclei, star 
formation is weaker, and their Wnp is beyond that at- 
tainable by CSF. The truncation time scale indicated by 
these measurements for these bins are from 80 to 100 
Myr. Galactic dust contributes ~ 0.16 magnitudes of 
extinction at the redshifted wavelength of Ha. 

IRAS20087-0308: The morphology of the K-band im- 
age shows a single nucleus, implying that the two galaxies 
have fully merged. Similarly we see one centrally located 
position on the slit with continuing star formation. The 
SNR for this object makes Whp difficult to measure, 
but the nuclear measurement in this object does appear 
slightly lower than the extended regions. This difference 
is not significant enough to measure a gradient in the 
population age, so we estimate an overall age of 100 Myr 
since star formation stopped across the disk. Galactic 
dust contributes 0.31 magnitudes of extinction at the 
redshifted wavelength of Ha. 

IRAS17208-0014: This galaxy has a single K band 
nucleus with many diffuse disturbed regions around it 
in the R band. The aperture with peak star formation 
(0.64 ± 0.03 Mq yr-i kpc'^) is adjacent to the region 
with the smallest Wh/s (4.7 ± 0.2 A) which corresponds 
to 170±30 Myr of continuous star formation (Figure[2le). 
Next to the K band centroid at the max of the R band 
continuum Whp — 6.6 ± 0.3 A and the star formation is 
diminished to 0.016 ± 0.003 Mq yr"! kpc'^. This aper- 
ture shows an underlying population with a truncation 
time scale of 57 ± 10 Myr. The furthest aperture with 
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a center 4.25 kpc from the nucleus has a shghtly larger 
and very little star formation, making the time 
since truncation of this population slightly larger at 100 
± 35 Myr. Galactic dust contributes ^ 0.67 magnitudes 
of extinction at the redshifted wavelength of Ha. 

IRAS18368+3549: The K-band morphology is a sin- 
gle nucleus with the R band revealing a faint tidal region 
encircling the nucleus. From this morphology we expect 
that the nuclei have completely merged and tidal distor- 
tions have had tirne to circularize. The nucleus has an 
Wh0 of 6.8 ± 0.4 A putting it beyond the maximum pos- 
sible in CSF. This measurement represents a time scale 
of 77 ± 17 Myr in the truncated constant star formation 
history, and 95 ± 17 Myr in the 6-in model. Whi3 in the 
extended region at 10 kpc from the nucleus is very large 
putting it beyond the model describable with these star 
formation histories. 

IRAS23365+3604: The morphology of this galaxy is 
similar to that of IRAS18368+3549, in that it has a 
bright nucleus with a tidal region that has had time to 
circularize. The nuclear regions have a peak in star for- 
mation rate (2.60 ±0.02 Mq yi'^ kpc"^), and low Wnp 

(3.3 ± 0.1 A). We measure 45 ± 5 Myr since continuous 
star formation began. Apertures exterior to this have 
similar times since a truncation of star formation (37 ± 3 
Myr in a truncated constant star formation history, 52±4 
Myr in a 6-in model). At a distance of 5.3 kpc from the 
nucleus, the slit intersects an extended tidal region, the 
center of which shows a dip in equivalent width (3.9 ±0.4 
A) in an otherwise increasing Whi3 with radius. In the 
furthest aperture no Ha emission is measured, and the 
measured Wnp of absorption indicates that truncation 
occurred 140 ± 25 Myr ago. Galactic dust contributes 
~ 0.19 magnitudes of extinction at the redshifted wave- 
length of Ha. 

IRAS09111-1007: This galaxy has undisturbed tidal 
regions characteristic of post-merger morphology (Figure 
[21a). The central kpc has Wh/} = 4.5 ± 0.2 A, which 
is indicative of younger stellar populations. From this 
point, increases toward the edge of the object, up 

to 9 -10 A at projected distances of 3 - 5 kpc. These 
large equivalent widths imply that 300 Myr has passed 
since the truncation of star formation. 

IRAS10378+1109: This galaxy has a compact nucleus, 
but it is measurable in several apertures. The fit of the 
absorption trough in the central aperture yields an un- 
physically narrow fit due to contamination by strong 
emission and is rejected from the rest of the analy- 
sis. This strong emission indicates strong star formation 
(6.73 ± 0.06 Mq yr-i kpc'^). The aperture 4.7 kpc has 
W/f;3 — 6.9 ± 0.4 Aindicating a truncation time scale of 
- 90 Myr. 

IRAS11506+1331: Star formation is maximum (2.23± 
0.06 M© yr~^ kpc~^) in the central most aperture (1.3 
kpc) . If a continuing star formation history were consid- 
ered, the absorption (W^^ = 6.2±0.25 A) would indicate 
that constant star formation has continued for 1 Gyr, 
however the presence of a strong absorption line indicates 
that there was a previous burst of star formation that 
dominates the continuum over the stars being produced 
in the current epoch of star formation. In this case, the 
end of the previous burst of star formation would have 



ended 69 ± 9 Myr ago. Exterior to this, regions have 
diminished star formation rates (< 0.1 M© yr~^ kpc~^) 
and truncation time scales up to 150 ± 35 Myr at 7 kpc. 

IRAS10565+2 U8: This object is classified as a mul- 
tiple merger by (jMurphv et al.l[r9 96V with a third com- 
ponent separated by 20 kpc from the brightest knot in 
our measurement. These data sample two of the nuclei 
and measure an increase in Wh/3 with radius. In the 
stronger nucleus Wnp drops to 1.8 ± 0.1 and has a large 
amount of star formation (7.28 ± 0.08 Mq yr~^ kpc~^) 
indicating a very young population. Wnp increases to 
8.3 ± 0.5 A, corresponding to a truncation time scale of 
^ 180 Myr. A dip in the trend occurs off center of the 
weaker nucleus. The truncation timescale in this region 
is ^ 50Myr. 

From analysis of the spatially-resolved sub-sample, 
we draw the following conclusions. Objects such 
as IRAS20087-0308, IRAS18368+3459, IRAS23365-3604 
and IRAS10565+2448 have regions at large radii that 
show no Ha, indicating that the star formation has been 
completely shut off. The others show examples of con- 
tinuing, but decreased, star formation activity at larger 
radii. Often the star formation rate per unit area drops 
by at least a factor of 10 at the largest radii measured. 
Though the gradient varies among the galaxies, we see 
regularity in that the time since truncation increases out- 
wards in radius. 

3.3. Other Sources of Spatial Gradients in Wnp 

We considered other physical effects that could cause 
large W^^^ values. Previous studies (iPoggian ti fc Wul 
l2000f ) attributed large Balmer equivalent width to in- 
termediate age stars having time to drift away from the 
dusty molecular clouds in which they were born. This 
effect can have a major effect in integrated spectra of en- 
tire objects. For dust to be the cause of larger equivalent 
width at larger radii, the population of hidden young 
stars must grow with distance from the nucleus. This 
conclusion contradicts measurements of the molecular 
gas distribution in ULIRGs. An active nucleus may con- 
tribute to the continuum emission in the central aper- 
ture. AGN typically contribute 15 — 30% of the total 
bolometric lum inosity from cool HII and LINER ULIRGs 
IVeilleux et al.l 12009). Boosting the continuum, however, 
only decreases the apparent stellar Whi3, so an AGN can- 
not explain large Whp values. Moreover, the gradients 
that we find on spatial scales larger than the seeing can- 
not be caused by an AGN. 

3.4. Comparison of Wnp among Nuclei 

Among the 25 spectra examined, the surface bright- 
ness of the merger remnant was not always high enough 
to measure the spatial dependance of the line profile. We 
can, however, compare Wnp among all the nuclei as il- 
lustrated in Figure [71 In the pre-merger subsample we 
measured Whp of both nuclei seen. Of these 18 nuclei, 4 
of them showed broad H/3 emission that did not allow the 
measurement of an underlying absorption line. A typical 
nucleus in the pre-merger sample presents Wh0 of 4 A, 
although the full range in W^^^ is large. In the merger 
sample, the median Whis is higher. Our sample is too 
small to be certain that this difference is significant; but, 
if verified, would suggest the central star formation de- 
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Fig. 6. — Three stellar population age measures are compared us - 
ing solar metallicity burst models fromlBruzual fc Charlotl II2003I ). 
Wji/^shows slightly larger equivalent widths than Wjy^. The 4000 
Angstrom break (D4000) is defined bv lKauffmann et al.l l|2003l ) as 
the ratio of fluxes in a bin redward of the break (4000 - 4100 A) to 
a bin blueward of the break (3850 - 3950 A). D4000 is monotonic 
and has an enhanced sensitivity over the range that the Balmer 
absorption function turns over creating an ambiguity in age. Mea- 
surements of D4000 as a function of position for these objects will 
more clearly define population age. 

clines once the nuclei merge. In the post merger phase, 
Whp would then continue to increase up to roughly 8 or 
9 A before starting to decline. 

4. DISCUSSION 

By resolving the recent star formation activity in 
ULIRGs over scales of 1-10 kpc, we learned that star for- 
mation is being turned off in ULIRGs from the outside 
inwards. Specifically, we measure an increase in Whi3 
with increasing radius and argue that the only reason- 
able interpretation is an increasing fraction of Type A 
stars (relative to earlier types). We go one step further 
here and interpret the reduced star formation rate as di- 
rect evidence for gas depletion. Above a threshold den- 
sity of approximately ~ 10 Mq pc~^, the star formation 
intensity empirically scales non-linearly with the gas sur- 
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Fig. 7. — Histograms of the nuclear Wh/B measurements are pre- 
sented for comparison between the different morphological classi- 
fications. The pre merger histogram does not include 4 of the 18 
nuclei that are present in the spectra due to the presence of emis- 
sion broad enough to make an absorption component can not be 
discerned. The merger stage objects are predominantly between 
4 and 7 A. The sample of old mergers is small, but the measured 
nuclei are greater than 5 A. 

face de nsity, 'Ssfb. /i^'^, where is the surface den sity 
of gas (jMartinfc KennicuttI 120011 : lKennicuttlll998[) . If 
star formation were the primary mechanism removing 
the gas, then the gas would be consumed more quickly 
in the highest density regions, r ^ fx/SFR cx 
i.e. the nuclei. Star formation alone produces a gradi- 
ent oppositely sloped compared to the one we measured. 
Our result therefore provides direct evidence that an- 
other mechanism moves gas out of the outer disk. It is 
natural to ask whether merger-induced infiows are that 
mechanism. 

Simulations of gas-rich, galaxy - galaxy mergers pre- 
dict an increase in star formation during the first pas- 
sage and a s econd, stronger starburst a.t the time of ac- 
tual merger CMihos fc Hernguistl [1991 ILotz et al.l 120081: 
Hopkins et al. 2009). The torques imparted on the gas 
during the merger cause the gas to flow inward, robbing 
the outer regions of fuel to generate stars. The infall is a 
well known result of the separation of stars from gas in 
the merger remnant; a consequence of the different coUi- 
sional properties of stars and gas (Hopkins et al. 2009). 
This separation allows the stars to impart a torque on 
the gas and cause infall. Models imply the correlation 
between the dynamical a ge of the merger a nd the age of 
the starburst population (|Lotz et al.ll200"8l) but have yet 
to predict how much of the star formation occurs in the 
central kpc. A key aspect of merger models should be 
to trace gas migration and the location of star formation 
activity througho ut the merger. 

Observations in iKewlev et al.l ()2006l ) provide circum- 
stantial evidence relating gas inflow to the starburst and 
therefore to the interaction. Galaxy mergers with greater 
starburst strength have lower nuclear metallicity. This 
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effect is expected to be due to the infall of pristine gas 
from the outer disk which dilutes the concentration of 
metals in the nucleus. We have measured the time of gas 
depletion at radii from 1 to 10 kpc. The results constrain 
the net infall velocity. 

4.1. Implication of Radial Gradients in Star Formation 

History for Gas Inflow Timescales 

Table [2] lists the time elapsed since the star forma- 
tion rate or, equivalently, the gas surface density dropped 
abruptly. The time elapsed since peak star formation ac- 
tivity is typically 100 to 300 Myr at 5 to 10 kpc falling to 
less than 50 Myr within 5 kpc. These timescales indicate 
the minimum time for gas inflow since the gas need only 
flow inward, not necessarily all the way to the center of 
the merger. 

In observations by iColina et al.l ()2005[ ). the 
dynamical mass for objects included in our 
study (IRAS17208-t-0014, IRAS20087-0308, IRAS 
23365-1-3604) are measured as similar to the Milky Way 
mass. The representation that we give for the orbital 
velocities may be reparameterized by substituting in 
distances and orbital velocities for other galaxies, as 
well as changing the density values used in the freefall 
calculation. Observations ('Dasvr a et al.l 120061 ) have 
indicated that ULIRGs are major mergers between 
galaxies with an average mass ratio of 1.5:1, indicating 
that the orbital velocities in each component are of 
the same order. This information allows us to choose 
the Milky Way mass value as the fiducial mass scale; 
and we estimate maximum infiow speeds of order 
v^nflow - 68 km s~\R/7 kpc)/(i*/100 Myr). 

We compare this to two timescales in ULIRGs: (1) 
the free fall timescale at a given radius in an isother- 
mal sphere and (2) the orbital period. Both of these 
timescales increase linearly with radius. To determine 
the freefall timescale we calculate the average density 
within particular radii using a density profile of the form 
Po {r/ro)-"^ with p„ = 1.3 x 10^ M p, kpc"^ the local 
halo mass density ()Gates et al.lll995l) . The average den- 
sity of the Milky Way interior to a radius of 8 kpc is 
Pmw « 4 X 10'' M0 kpc~^ giving the free fall timescale 
Tff{r) ~ 40 Myr PMw/p{r)- The orbital period is 
TorbiR) ~ 220 Myr (i?/8 kpc)/(w/220 km s'^). In Figure 
m we show the truncation timescale as a function of dis- 
tance from the center of the merger. The scatter in val- 
ues does not distinguish between freefall and orbital time 
scales. They do rule out the much longer time scales as- 
sociated with diffusion and lend em pirical support to the 
timescales assumed in simulations ()Hopkins et al.l [20081 
[2009h . 

4.2. Picture Emerging from Comparison of Stellar and 

Dynamical Ages 

Gas migration toward the center slows down star for- 
mation in the outer radii. The paucity of massive stars al- 
lows Whi3 to grow. Regions closer to the nucleus are fed 
infalling gas from exterior regions, prolonging star for- 
mation and allowing Whis to remain low longer. As the 
merger advances, the cessation of star formation trails in- 
ward. All 25 ULIRGs have star formation in their center, 
indicating the gas surface density there remains above 
the threshold density. As discussed in Section 13. 4[ the 
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Fig. 8. — The grey points are the measured values for age and 
position in the truncated star formation history model. The black 
points are the median values of age within bins that contain equal 
numbers of points, where the bin size is represented by the error 
bar. The median points show that the infall timescale is on the 
same order as the orbital period and the free fall time scale. The 
scatter at large radii may indicate differences in star formation 
history or infiow time scale. 

higher nuclear Wh/i in the merger and post-merger ob- 
jects may indicate the central star formation rate is de- 
clining by this stage. This result provides observational 
evidence that the central gas density starts declining once 
the nuclei have merged. 

5. CONCLUSION 

The physical processes regulating the rate of gas infall 
during mergers is important because it influences the age 
and metallicity gradients in the merger remnant as well 
as the strength of supernova and AGN feedback. By ex- 
amining the positional dependence of stellar spectral in- 
dices in ESI long slit spectra , we determined the recent 
star formation history across galaxy - galaxy mergers. 
Strong II/3 absorption indicates a diminished star forma- 
tion rate over the past few hundred Myr, which should 
be accompanied by an increase in D4000, while hydrogen 
Balmer emission indicates the presence of massive stars 
from more recent star formation. 

We find the measured Wnp increases with radius in 
the sample of spatially resolved objects. We attribute the 
large WHfs to a rapid decrease in star formation activity 
over 100 Myr ago at radii greater than 5 kpc. At radii 
of a few kpc, the activity appears to have decreased just 
50 to 100 Myr ago based on the slightly lower, but still 
prominent, Wh/s- Our measurements in the central kpc 
are consistent with star formation continuing unabated. 
The nature of the gradient implies that gas was removed 
from larger radii first allowing the stellar populations to 
age. 

We propose that these age gradients trace the inward 
flow of gas from large radii, adding to the central gas sup- 
ply and prolonging the central starburst. Our data pri- 
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marily address the pre-merger and early merger phases. 
At later stages, the smfaee brightness of a ULIRG is 
typically too low to make the equivalent measurement. 
We do, however, find preliminary evidence that these 
spectral indices in the nuclei of late mergers indicate de- 
creasing star formation. We interpret these results as 
direct evidence for strong gas inflow in the pre-merger 
and early merger phases with a possible suppression by 
the late merger phase. 

Assuming the time since star formation diminished re- 
flects the gas inflow timescale, our measurements indi- 
rectly constrain the gas infall rate. The time scales for in- 
fall to occur are within the range of the orbital timescale 
and the freefall timescale of an isothermal sphere. 
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TABLE 2 
Position Dependent Measurements 



IRAS name v — fnuc 

(1) 



(2) 



(3) 



CSF 



(4) 



TSFH 



(5) 



DSFH 



(6) 



SFR 



-7.0 ± 0.8 


8.8 


it 


0.5 


N/A 


-5.0 ± 0.8 


7.0 


± 


0.3 


N/A 


-1.9 ± 1.0 


4.8 


± 


0.1 


203.4 


0.5 ± 0.7 


5.0 


it 


0.1 


221.3 


2.2 ± 0.7 


6.1 


± 


0.2 


895.3 


4.6 ± 0.8 


6.2 


it 


0.2 


N/A 


7.3 ± 1.0 


8.4 


it 


0.6 


N/A 


-8.6 ± 0.9 


7.0 


it 


0.5 


N/A 


-6.4 ± 0.8 


4.7 


it 


0.2 


182.0 


-4.5 ± 0.7 


4.7 


it 


0.2 


182.7 


-2.7 ± 0.7 


4.9 


± 


0.2 


209.6 


-0.8 ± 0.6 


5.6 


it 


0.2 


360.0 


0.9 ± 0.7 


7.0 


it 


0.2 


N/A 


3.1 ± 0.9 


6.6 


it 


0.4 


N/A 


-4.9 ± 1.0 


8.5 


± 


0.9 


N/A 


-2.4 ± 1.0 


7.3 


± 


0.7 


N/A 


-0.4 ± 0.6 


7.3 


± 


0.8 


N/A 


1.9 ± 1.0 


6.3 


it 


0.8 


N/A 


4.6 ± 1.2 


6.7 


it 


1.2 


N/A 


-4.2 it 0.4 


7.2 


it 


0.6 


AT / A 
JN / A 


Q 1 -1- n /I 
-o.i it U.4 


6.6 


± 


0.3 


IN / A 


-2.1 ± 0.4 


6.4 


± 


0.3 


N/A 


-1.1 ± 0.3 


5.2 


it 


0.4 


260.1 


-0.3 ± 0.3 


4.7 


it 


0.2 


174.0 


0.6 ± 0.3 


6.2 


it 


0.4 


N/A 


1.8 ± 0.4 


12.0 it 0.9 


N/A 


-10.0 ± 1.8 


12.9 it 1.6 


N/A 


-2.5 ± 0.9 


9.6 


it 


0.7 


N/A 


-0.4 ± 0.7 


6.8 


± 


0.3 


N/A 


2.4 it 1.3 


6.8 


± 


0.4 


N/A 


11.4 it 2.1 


7.5 


it 


1.4 


N/A 



152454-1019 



20046-0623 



20087-0308 



17208-0014 



18368-t3549 



-1-13.8 
-20.7 
-H23.7 

-15.5 
+404.4 

-404.4 



-23.8 
+33.0 
-33.1 
+63.2 
-73.0 



+71.2 
-71.2 
+21.0 
-14.3 
+ 1.6 
-1.4 
+ 1.9 
-1.8 
+8.4 
-8.3 
+8.1 



250.1 

88.5 

23.3 

25.4 

48.6 

53.1 



88.1 
21.8 
21.9 
24.0 
34.6 
88.5 
67.1 



+34.5 

-26.2 

+3.1 

-3.0 

+2.7 

-2.6 

+3.2 

-2.7 

+3.8 

-3.7 

+ 16.1 

-10.5 

+ 19.6 

-16.8 



187.i 

106.7 

107.1 



+65.1 
-65.1 
+46.7 
-42.0 
+49.4 
-44.0 

55.7 tf,\ 
73.9 

100.4 ^36.0 

+ 13.7 
-12.9 
+ 11.7 
-9.3 
+5.9 
-5.4 
+2.5 
-2.5 
+ 12.6 
-11.0 



70.2 
56.8 
28.4 
20.9 
49.8 
N/A 



N/A 
N/A 
76.7 
79.6 
123.8 



+ 17.1 
-16.3 

+26.9 
-21.2 



^68. 3 
-68.3 
+23.8 
-14.0 
+ 1.1 
-0.7 
+ 1.3 
-1.4 
+8.3 
-7.2 
+8.3 
-7.9 

197.8 tfi^ 



263.6 
106.8 
38.i 
40.3 
66.5 
70.9 



+37.5 
-26.8 
+ 1.9 
-2.7 
+ 1.6 
-2.4 
+2.3 
-1.6 
+5.3 
-5.3 
+ 16.7 
-10.4 
+19.8 
-16.2 

+60.0 
-60.0 

127 7 +*''''^ 
-45.4 

128.2 1^6.9 
7Q c +36.5 

I 'J•^J 22 \ 

qn r +70:9 
y^-O _42.1 



106.5 

38.0 

38.1 

39.2 

50.4 

106.8 

85.3 

204.1 



116.6 tfd 
+ 12.9 

-13.9 
+ 12.4 
-9.2 
+ 7.8 
-3.7 
+ 1.6 
-2.2 
+ 12.4 
-10.8 



89.0 
74.6 

42.4 
37.2 
67.6 
N/A 



N/A 
N/A 

qr 9 +15.3 
WJ.^ -16.9 

9^.9 _2i.6 

145.11^1^1^ 



0.0014 it 0.0007=^ 
0.026 it 0.016 « 
1.52 it 0.02 
4.83 it 0.05 
1.01 it 0.05 
6.8 it 2.1 
0.008 it 0.006 * 

0.0053 it 0.0005 
0.156 it 0.002 
3.95 it 0.05 
1.18 it 0.01 
3.62 it 0.08 
8.3 it 0.7 
0.98 it 0.32 



0.058 it 0.025^ 

d 

0.8 it 0.2 
0.04 it 0.05 

0.004 it 0.003^ 
0.016 it 0.003 
0.099 it 0.004 
0.52 it 0.02 
0.51 it 0.01 
0.64 it 0.03 
0.11 it 0.02 



0.8 it 0.5 ^ 

d 

0.76 it 0.26 
0.00008 it 0.00009 ■ 



23365+3604 



-8.7 it 0.8 


7.8 


it 


0.4 


-6.8 it 0.5 


7.0 


it 


0.4 


-5.3 it 0.6 


3.9 


it 


0.4 


-3.7 it 0.7 


5.6 


± 


0.3 


-1.9 it 0.6 


4.9 


it 


0.2 


-0.5 it 0.4 


3.4 


it 


0.1 


0.7 it 0.5 


3.3 


it 


0.1 


1.7 it 0.4 


5.7 


it 


0.1 


2.8 ± 0.4 


6.7 


± 


0.2 


4.0 it 0.4 


5.3 


it 


0.2 


5.4 it 0.7 


5.9 


it 


0.6 



N/A 
N/A 

07 7 +42.9 
-29.0 

358.4 

209.0 tjs ? 

50.0 l^:^' 

46.4 
392.1 

N/A 
277.0 
N/A 



141.2 
87.5 
11.4 
34.5 
23.9 
2.5 
2.7' 
36.9 
75.8 
30.0 



+ 29.6 
-24.6 
+27.3 
-19.8 
+4.9 
-6.2 
+4.8 
-4.5 
+2.6 
-2.2 
+3.0 
-1.0 



+2.5 
-2.5 
f9.2 
-10.4 
+3.4 

40.7 ig^ 



^26. 8 
-23.4 
+30.1 
-19.9 
+3.4 
-25.7 
+ 7.4 
-6.3 
+ 1.9 
-1.2 

4.1 +°i 



161.2 

105.9 

29.9 

50.3 

39.1 



-0.2 



2.7<= 
52.3 
94.3 
44.0 
60.0 



+5.7 
-2.1 

+8.8 

-11.1 

+4.7 

-2.8 

+19.0 

-14.8 



b 
b 

0.011 it 0.001 
0.036 it 0.011 
0.20 it 0.01 
1.39 it 0.02 
2.60 it 0.02 
1.83 it 0.06 
1.27 it 0.37 
0.025 it 0.023^* 



09111-1007 



-3.2 it 0.6 


8.8 it 0.9 


N/A 


-1.9 it 0.4 


7.4 it 0.4 


N/A 


-0.4 it 0.4 


6.5 it 0.2 


N/A 


1.1 ± 0.4 


4.5 ± 0.2 


150.8 


2.1 it 0.3 


6.3 it 0.2 


N/A 


3.1 it 0.4 


7.6 it 0.3 


N/A 


4.2 it 0.4 


7.5 it 0.5 


N/A 


5.5 it 0.5 


10.0 it 0.8 


N/A 


7.0 it 0.6 


10.0 it 1 


N/A 


-6.7 it 1.9 


7.7 it 1.2 


N/A 



+21.8 

-18.0 



243 fi +95-8 
zy:O.D _95 8 

+ 21.5 

-24.4 

+ 11.2 

-8.2 

+2.2 

-2.1 

+8.0 

-7.5 

+ 17.1 

-16.3 

121.1 t'^l 

N/A 

N/A 

133.0 



257.6 



+90.4 



112.4 

64.1 

18.6 

55.2 

128.5 



-90.4 

133.5 tfA 

81.6 t^;^ 
+2.0 
-1.6 
+7.6 
-7.4 

+ 15.7 
-16.1 

142.5 

N/A 
N/A 

153 7 +121.6 
100.1 _74,9 



35.1 
73.0 
149.5 



0.19 it 0.04 
0.90 it 0.07 
1.35 it 0.03 
0.171 it 0.005 
0.049 it 0.001 
0.001 it 0.0001 
0.0003 it 0.0001 
0.00008 it 0.00005'' 
0.00001 it 0.00001'' 

0.0073 it 0.0025 
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TABLE 2 



Continued 



IRAS name 


(1) 


(2) 


CSF 

(3) 


TSFH 

(4) 


DSFH 

(5) 


SFR 

(6) 




-1.7 ± 1.4 

1.3 ± 1.3 

4.4 ± 1.0 
7.1 ± 1.1 
10.4 ± 1.3 


C 

6.2 ± 0.2 
7.1 ± 0.4 
7.7 ± 0.5 
6.7 ± 0.4 


N/A 
N/A 
N/A 
N/A 
N/A 


N/A 
51.4 tli 
94 6 

131 1 +35.5 
79 o +21.1 


N/A 
69.2 +|-? 

Ill 1 +?2 n 
^^^•^ —18.9 

151.9 lii 

-'l-'J -19.1 


0.60 ± 0.01 
2.23 ± 0.06 
0.053 ± 0.007 
0.015 ± 0.001 
0.028 ± 0.001 


10378+1109 


-4.4 ± 1.7 
-0.2 ± 1.1 
2.4 ± 0.9 
4.7 ± 1.1 


3.6 ± 0.4 

e 

3.9 ± 0.2 
6.9 ± 0.4 


62 6 +38.7 

DZ.D _27.8 

N/A 

Q1 Q +16.0 

-14.6 

N/A 


fi Q +6. .5 
"•"^ -12.4 

N/A 
11.9 11° 

OO.U _20.8 


A Q +26.3 
^•■^ -0.4 

N/A 
30.2 11;° 
103.2 tliA 


2.5 ± 1.2 
6.73 ± 0.06 
1.04 ± 0.05 
0.37 ± 0.20 



10565-1-2448 



-7.5 ± 0.6 
-6.3 ± 0.5 
-4.9 ± 0.5 
-3.4 ± 0.5 
-2.3 ± 0.5 
-1.2 ± 0.4 
-0.1 ± 0.5 
1.2 ± 0.6 
2.2 ± 0.4 



7.2 ± 0.2 
8.1 ± 0.2 
6.1 ± 0.6 

8.3 ± 0.5 
6.6 ± 0.3 
4.3 ± 0.2 
1.8 ± 0.1 
5.8 ± 0.2 
7.8 ± 0.4 



N/A 
N/A 

884.4 tint 
N/A 
N/A 
132.2 
N/A 
425.7 tll'e' 
N/A 



+ 24.0 
-22.6 



+ 14.7 
-15.3 
+ 11.4 
-11.8 
+24.3 
-14.2 
+72.4 
-38.6 
+ 16.8 
-13.1 
+2.6 
2.3 



100.1 
158.3 
48.5 
176.0 
65.7 ; 
16.5 
N/A 
38.5 tlr 
140.9 tin 



66.5 



1 ic 1 +19.9 
116.1 1 

1 7fi 7 +10.3 

I'"-' -10.6 

+25.2 

-16.3 

1Q9 7 +69-5 
LUZ.I _35 

+17.0 

-13.2 
+2.1 
-2.1 



83.6 
33.5 
2.7<^ 
55.9 



+6.6 

160.9~li:° 



b 
b 
b 

0.010 ± 0.008^ 
0.34 ± 0.13 
0.42 ± 0.01 
7.28 ± 0.08 
2.8 ± 0.1 
0.14 ± 0.08 



Note. — Col.(l): Projected distance from nucleus to center of measured region in kiloparsccs. The error corresponds to the 
spatial width of the aperture. Col. (2): Absorption equivalent width for the aperture measured in A . Col. (3): Stellar population 
age in continuing star formation history measured in Myr. In the apertures where the Wjf/j is too large for the given stellar 
population, the section is marked with "N/A". Col.(4): Stellar population age in Truncated star formation history measured 
in Myr. Col. (5): Stellar population age in S function star formation history measured in Myr. Col. (6): Star formation rate 
per aperture measured in M© yr"^ kpc~^. 

In these apertures, H/3 emission was not detected, indicating heavy extinction. These values represent a lower limit to the 
SFR based on assuming that the H/3 emission line has approximately the same peak as the noise, making it unmeasurable. 
^ Ha emission was not detected in these apertures. 

Upper limit of age is presented for very low Whb- 

These apertures have forbidden/Balmer line ratios that indicate a considerable contribution by AGN. Star formation rate is 
excluded due to possible confusion. 
The fit provided an unphysically narrow fit of absorption due to contamination from emission, so was rejected from analysis. 



